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We present results of numerical simulations for pure U(\) gauge theory in a non-commutative 
space. The theory is mapped onto a dimensionally reduced matrix model, which renders its nu- 
merical treatment feasible. New data on large lattices reveal the scaling of Wilson loops and their 
correlation functions in the simultaneous limit to the continuum and to infinite volume, at fixed 
non-commutativity. In this on-going project we are particularly interested in the IR behaviour, the 
"photo-ball" spectrum and in the photon dispersion relation. 
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1. Non-commutative U(l) gauge theory 

Somewhat overshadowed by the celebrations for Einstein's works of 1905, we are also cele- 
brating this year the 200th anniversary of Sir William Rowan Hamilton [jj]], particularly here in Ire- 
land. One of his achievements was the discovery of quaternions in 1843: he had thought for a long 
time about possibilities to extend the representation of complex numbers by two real components 
to a system with three real components, until he noticed that he had to proceed to four components 
to arrive at a sensible system — the first non-commutative (NC) algebra that was studied. 

In the 20th century, the concept of non-commutativity for space coordinates and momenta 
became standard in quantum physics, but also the idea of NC space coordinates is about 60 years 
old [g]. It experienced a powerful renaissance in 1998, triggered by the identification of open 
strings at low energy with NC field theory [[|]. NC spaces are also considered a promising approach 
to quantum gravity [QJ. On the phenomenological side, it might for instance explain the observation 
of high energy photons from far away galaxies, beyond the energy threshold apparently predicted 
by the Standard Model Here we study NC gauge theory in its own right. 

In that framework, quantum mechanical position operators obey a commutation relation of the 
form [jCju, jc v ] = i®nv , where we assume the non-commutativity tensor & to be constant in (Eu- 
clidean) space-time. More precisely, we consider the case of two commutative directions (which 
include the Euclidean time), and an NC plane with the relation [xi,Xj] = idEjj (ij 6 {1,2}). 

Field theory on such a space can be written in terms of our usual (commutative) coordinates x, 
if all field multiplications are performed by star products, 

(j)(x)*\i/(x) := 0(x) exp (-^06^ ~^/) . (1.1) 
In particular, the action of pure U(l) gauge theory takes the form 

S\A] = Jd 4 x Tr[F MV *F MV ] , = d^A v - d v A^ + i[A^*A v -A V *A M ] , (1.2) 

which is star-gauge invariant. The Yang-Mills type self-interaction term is expected to yield a 
"photo-ball" spectrum [Bj .It may modify the photon dispersion relation at low energy, as it was 
observed non-perturbatively in the NC A0 4 model [Q], as a consequence of UV/IR mixing effects 
[§]. These effects drastically complicate the perturbative treatment. However, the j3 function could 



be computed, suggesting asymptotic freedom [gp. There are a number of further perturbative [ |10| ] 



and semi-classical [11] studies 



A formulation on a (fuzzy) lattice is possible [J12J], and it relates the spacing a on a. N x N 
lattice to the NC parameter as 

d = -Na 2 . (1.3) 

71 

Then the Double Scaling Limit (DSL), which takes simultaneously N — > °° and a — > at 6 = 
const. , leads to a continuous NC plane of infinite extent. 



'The term "photo-ball" is an obvious analogue to the glueball; to our knowledge, it first occurred in the work by 
Fatollahi and Jafari quoted in Ref. W. 
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But such a formulation is not immediately applicable for simulations, in particular because of 
the request for star-unitary link variables. However, there is an exact map Jl3| ] of the NC N x N 
lattice onto the dimensionally reduced Twisted Eguchi-Kawai model [14] with the action 



S[U] = -fiN £ %j Tr [UiUjUjuJ] , 3T U = 3% x = exp{ni{N + l)/N} , N odd , (1 .4) 

¥J 

where £/; € U (N), i = 1,2. The (analogue of a) rectangular Wilson loop of sides al and a J is now 
given by 



W ij (IxJ) = ^ffTr(ufujU? I ul J ) . 



(1.5) 



Mapping this quantity back to the lattice yields in fact a sensible definition of a Wilson loop in the 
NC gauge theory [JlS]|. Note that NC Wilson loops are star-gauge invariant and complex. 



2. Numerical results for the Double Scaling behaviour 

This mapping of the NC plane onto a matrix model (one in each lattice site of the commutative 
plane) enables numerical simulations. The next challenge is to identify the dimensional lattice 
spacing a(/3) in order to evaluate observables in the DSL. The simple ansatz 

aoc\/p (2.1) 
turned out to be successful, as we are going to illustrate in a sequence of plots. Along with relation 



< JTT3| ) it implies N/j5 2 = const. 2 We always deal with N 2 x(N± l) 2 lattices, where N 2 ( (N±\) 2 ) 
is the lattice size in the NC (in the commutative) plane, and the NC plane is mapped onto a twisted 
Eguchi-Kawai model. We present results at N = 45, 55, 65, 71 and 81, and the corresponding j3 
values are fixed such that N/f} 2 = 20 in all cases. 3 This means that we are always in the weak 



coupling phase [14, 16] 



Figure [j] (on the left) shows the real part of the Wilson loop in the NC plane as a function of 
the loop area. We see that the loops of the same area I 2 a 2 reveal a convincing Double Scaling, if 



we insert the ansatz (|2JJ) (we chose the proportionality constant = 1). The same is true for the 



phase of these Wilson loops, as the plot on the right-hand-side of Figure [j] shows. Qualitatively this 



behaviour is similar to the Wilson loops in 2d NC QED Q17Q : small loops are almost real and decay 
exponentially as the area increases. On the other hand, for large loops the real part oscillates around 
zero, and the phase grows linearly. The latter property is reminiscent of the Aharonov-Bohm effect, 
if one identifies 6 with an inverse magnetic field across the NC plane, as suggested by Peierls Jig]. 

Next we consider correlation functions of Wilson loops, which are located in the NC plane but 
separated in the (commutative) Euclidean time. As examples, we show the correlation of the real 
parts of 4 x 4 loops (on the left), and of the phases of 2 x 2 loops (on the right). In both cases we 
still observe a decent Double Scaling. Its quality improves as N in increases, which indicates that 



our simple ansatz (2. 1 ) does indeed describe the DSL limit asymptotically. 



In contrast, in d = 2 we had to set a <>= 1//3 2 , and therefore N/fi = const. [0]. 

3 We are also working on a systematic search for Double Scaling by matching the data without any assumption about 
the relation between /3 and a. This represents a completely unbiased test of the DSL postulated here. 
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Figure 1: The Wilson loop in the NC plane: its real part (on the left) decays at small area, and for larger 
areas it oscillates around zero. In that regime, the complex phase (on the right) turns sizeable and it begins 
to grow in an (approximately) linear way. 
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Figure 2: The correlation of Wilson loops in the NC plane, separated by a distance aAt in Euclidean time. 
Our examples are the correlation the real parts of 4 x 4 Wilson loops (on the left), and of the phases of 2 x 2 
loops (on the right). 

Finally we also consider Wilson loops in the commutative plane. These loops, as well as those 
in the mixed planes, are real due to the reflection symmetry on the commutative axes. From Figure 
|3| (on the left) we see that for loops in this plane the convergence towards the DSL is a little more 
laborious, but for our largest ./V values it sets in also here. On the right-hand-side of this Figure 
we also consider the correlator of 6 x 6 loops in the commutative plane, again separated in the 
Euclidean time, which confirms the above statement with respect to the DSL. 



3. Conclusions and outlook 

We reported on our progress in a numerical investigation of pure (7(1) gauge theory in a NC 
space-time. In particular, we presented results which reveal the asymptotic rule for the Double 
Scaling, which takes the system to the continuum and to the infinite volume at the same time; 
the entanglement of these limits can be viewed in the light of the notorious UV/IR mixing. We 
observed that the simple ansatz ( |Q| ) for the lattice spacing works in a satisfactory way. 

Hence the bases is now provided to study the observables of physical interest and extrapolate 
them to the DSL. In particular, the Wilson loop correlators — examples are shown in Figure || 
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Figure 3: The Wilson loop in the commutative plane: the plot on the left depicts the loops themselves (they 
are real) against the dimensional area. On the right we show the correlation of 6 x 6 loops, separated in 
Euclidean time. 

— provide the bases for the evaluation of the "photo-ball spectrum". We also hope to obtain 
results for the NC distortion of the photon dispersion relation, which could then be confronted 
with experimental data. Several high precision experiments dealing with cosmic rays are about 
to measure the photon dispersion to a very high accuracy, over a broad range of energies (see for 
instance Ref. [0]). 

While our results suggest the existence of a finite continuum theory, we note that perturbative 
calculations revealed an infrared instability of the trivial vacuum p0| ] (as long as the model is not 
rendered supersymmetric). Indeed, we do observe numerically that the open Polyakov lines (which 
are star-gauge invariant as well) acquire non-zero expectation values, in accordance with perturba- 
tion theory. We therefore consider that we are actually probing a stable vacuum, which might be 
obtained after the condensation of "tachyons" in the trivial vacuum. This issue shall be discussed 
in a forth-coming paper. 
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